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It has been observed experimentally [H.R. Xia, C.Y. Ye, and S.Y. Zhu, Phys. Rev. Lett. 77, 1032 
(1996)] that quantum interference between two molecular transitions can lead to a suppression or 
enhancement of spontaneous emission. This is manifested in the fluorescent intensity as a function of 
the detuning of the driving field from the two-photon resonance condition. Here we present a theory 
which explains the observed variation of the number of peaks with the mutual polarization of the 
molecular transition dipole moments. Using master equation techniques we calculate analytically 
as well as numerically the steady-state fluorescence, and find that the number of peaks depends on 
the excitation process. If the molecule is driven to the upper levels by a two-photon process, the 
fluorescent intensity consists of two peaks regardless of the mutual polarization of the transition 
dipole moments. If the excitation process is composed of both a two-step one-photon process and a 
one-step, two-photon process, then there are two peaks on transitions with parallel dipole moments 
and three peaks on transitions with antiparallel dipole moments. This latter case is in excellent 
agreement with the experiment. 



I. INTRODUCTION 



There have been a large number of theoretical studies on the effects of quantum interference in atomic and molecular 
systems This phenomenon was first suggested by Agarwal || who showed that the spontaneous emission from 
a degenerate V-type three-level atom is sensitive to the mutual orientation of the atomic dipole moments. If they 
are parallel a suppression of spontaneous emission can appear and a part of the population can be trapped in the 
excited levels. Similar predictions were reported for other configurations of three- and multi-level atoms and show 
that quantum interference can lead to many interesting effects such as amplification without population inversion || , 
electromagnetically-induced transparency [Q , phase dependent spectra and population inversions Q , and ultranarrow 
spectral lines ||. 

Zhu and Scully Q and Lee et al. || have shown that in the case of a non-degenerate V-type atom driven from an 
auxiliary level, quantum interference can lead to the elimination of the central line in the fluorescence spectrum when 
the driving field is tuned to the middle of the upper levels splitting. This interesting effect suggests that quantum 
interference can be used as a mechanism for controlling and even for suppression of spontaneous emission. 

In 1996, Xia et al. (|] carried out the first experimental investigation of constructive and destructive interference 
effects in spontaneous emission. In the experiment they used sodium dimers, which can be modeled as five-level 
molecular systems with a single ground level, two intermediate and two upper levels, driven by a two-photon process 
from the ground level to the upper doublet. By monitoring the fluorescence from the upper levels they observed 
that the total fluorescent intensity, as a function of two-photon detuning, is composed of two peaks on transitions 
with parallel and three peaks on transitions with antiparallel dipole moments. The observed variation of the number 
of peaks with the mutual polarization of the dipole moments gives compelling evidence for quantum interference in 
spontaneous emission. 

It is our purpose in this paper to present a theoretical explanation of the observed fluorescent intensity and, in 
particular, to explain the variation of the number of the observed peaks with the mutual polarization of the molecular 
dipole moments. We point out here that the previous theoretical studies [^J^] of quantum interference between two 
transitions with parallel or antiparallel dipole moments have dealt with the fluorescence spectrum. By contrast in 
the experiment, the total fluorescent intensity, as a function of two-photon detuning, was observed. Agarwal [Tc| ] 
has provided an intuitive picture for the observed spontaneous emission cancellation in terms of interference path- 
ways involving a two-photon absorption process. Recently, Berman [ [Ti| has shown that the experimentally observed 
cancellation of spontaneous emission involving a two-photon absorption process can be interpreted in terms of pop- 
ulation trapping. Although a cancellation of spontaneous emission is present with a two-photon excitation process, 
no variation of the number of peaks with the polarization of the dipole moments exist in the fluorescent intensity. 
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In summary, no explanation has been offered until now for the observed variation of the number of peaks in the 
fluorescent intensity with the mutual polarization of the transition dipole moments. 

In this paper we consider a five-level system driven by a single-mode coherent laser field, which models the ex- 
perimental configuration set up by Xia et al. B. Working with the master equation of the system, we calculate the 
steady-state fluorescent intensity as a function of the laser frequency for two different transitions from the upper levels 
to intermediate levels. One transition is in the visible region and has parallel dipole moments. The other transition is 
in the ultraviolet and has antiparallel dipole moments. We assume that there is spontaneous emission from the upper 
to the intermediate levels and thence to the ground level so the dynamics of the system are restricted to these five 
levels. In a real sodium molecule, the situation is more complex, with other decay channels, and laser-field couplings 
between various real states However, we believe that our simple model does explain the basic physical effects 

which have been observed in the experiment. 

In Xia's paper 0J the excitation of the upper states is described as a two-photon process. As we will see later, the 
two photon excitation process can only ever lead to two peaks in the fluorescent intensity, independent of the mutual 
polarization of the dipole moments. We show that the experimentally observed variation of the number of peaks arises 
from the presence of an additional two-step, one-photon excitation processes. 

The paper is organized as follows. The master equation for the five-level molecular system driven by a single mode 
laser field is derived and analyzed in section II. The analytical and numerical results for the total fluorescent intensity 
for the two-photon coupling only are studied in section III. In section IV, we investigate the corresponding results 
when the system has both one- and two-photon coupling. We also examine the approximations made and make 
comparisons with the experimental results. A discussion is given in the concluding section V. 



II. THE MASTER EQUATION 

The energy-level scheme of the system we are considering is shown in Fig. |l|, in which we follow the notation 
of Ref. !|. The five-level molecule consists of two upper levels |oi) and I02), two intermediate levels |6) and \d), 
and a single ground level |c). The upper levels are separated by the frequency w 12 which is much smaller than the 
frequencies uin and u>2b of the \a\) — ► \b) and |aa) — ► \b) transitions and the frequencies U)\d and u>2d of the \ai) — ► \d) 
and 1 02) — » \d) transitions. As in the sodium dimers used in the Xia's experiment 0], we assume that the frequencies 
u>ib and u>2b are significantly different from the frequencies uji^ and u>2d- The transitions |ai), I02) — > \b) correspond 
to the visible region, whereas the transitions \ai), {0,2} — > \d) correspond to the uv region. 

In the molecule, the one-photon transitions |cti), \a%) — > \b), \d) — > |c) are connected by electric dipole moments, 
whereas the transition \a±) — > \a2) and the two-photon transitions |cii), \a%) — > |c) are forbidden in the electric dipole 
approximation. The molecular dipole moments can have different orientations (polarizations) and two dipole moments 
which are close in frequency can interfere with each other if they are not orthogonal. In the experiment, a destructive 
interference was observed between two transitions, \a±) — ► \b) and |o2) — > \b), with parallel dipole moments, and a 
constructive interference was observed between transitions \ai) — * \d) and (02) ~~ * \d) with antiparallel dipole moments. 

In order to quantify the mutual orientations of the transition dipole moments, we introduce a parameter 

\HijWki\ 

where flij is the matrix element of the transition dipole moment between \i) and |j) levels. Using the subscripts u 
and v to denote the ultraviolet and visible transitions in the experiment, we have p u — 1 (parallel dipole moments), 
while p v = — 1 (antiparallel dipole moments). 

For simplicity we will assume that the magnitude of the interfering dipole moments are the same. Thus, the upper 
doublet decays to level \b) at rate j v = 71k = 726 and to level \d) at rate j u = 71^ = 72^. Here u and v again refer to 
visible and ultraviolet. The intermediate levels \b) and \d) decay to the ground level |c) at rates "fb and 7^ respectively. 

The system is driven by a single-mode tunable laser of frequency ojl- In the experiment the dye laser was coupled 
to the two-photon transition |c) — ► |ai), (02) m order to avoid the Doppler effect (which we ignore in our analysis). 
Here, we must ask a question whether the two-photon coupling in the experiment was the only coupling of the 
laser to the system. It is stated in the experimental paper 0J that the two-photon transition in sodium dimers was 
enhanced by a near-resonant intermediate level, indicating that the laser could also couple the ground state |c) to 
the upper states |ai), \a%) via cascaded one-photon transitions. Here, to avoid introducing an extra level, we take the 
near-resonant intermediate level to be \b), so the laser can also produce a two-step one-photon transition |c) — > \b) 
then I b) — > \a\), \a2)- In our opinion this channel of the excitation was possible in the experiment as the one-photon 
transitions in the molecule are in the visible region and their dipole moments are parallel ||. We will see later that 
the presence of this channel of excitation will be crucial in the explanation of the experimentally observed fluorescent 
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intensity profile. With only two-photon excitation quantum interference can be observed but the fluorescent intensity 
exhibits two peaks (as a function of laser detuning) regardless of the mutual orientation of the transition dipole 
moments. The three peak structure of the fluorescent intensity observed in the experiment in the uv region can only 
result from the presence of the two-step one-photon channel. 

We calculate the steady-state intensity of the fluorescence from the upper doublet to the intermediate levels as 
follows. The intensity is proportional to the normally ordered first-order correlation function of the scattered field 



I (f, t) cx (M^ (r, t) ■ E {+) (r, £)\ , 



(2.2) 



where (r, t) is the positive frequency part of the electric field operator at a point r in the far-field zone of the 
system outside the driving laser field. In terms of the density matrix elements of the system the scaled steady-state 
(t — ► oo) intensity on the ultraviolet and visible transitions is 



(pn + P22 + 2p u/v Rep 12 ) 



(2.3) 



Here pn and p 22 are the steady-state populations of the level |di) and {0,2), and pi 2 is the stead-state coherence 
between them. 

We find steady-state values of the populations and coherences from the master equation of the system. The master 
equation can be written in the Lindblad form [|13| as 



P = CrcvP + ArrP, 

where the reversible and irreversible terms are, respectively 

CrcvP = -i[H, p], 

and 

Arr = 7, (1 + Pv)T> [|6>((oi| + (oaD/VS] +7, (1 - Pv)V [|6)((oi| - (o 2 |)/V2 

+ 7« (1+Pu)© [[|d>«Oi| + <a 2 |)/V2] + lu (l-p u )v[\d){( ai \ - (a 2 \)/V2 
+ lb V[\c)(b\]+ ld V[\c)(d\] 
= lv V \\b) «oi| - (oa I )] + luV \\b) ((ai| + (a 2 1 )] + lb V[\c) (6|] + ld V[\c) (d\] . 

Here, T> is a superoperator defined for arbitrary operators A and B as 

V[A]B = ABA^ - i{AU, B}. 



(2.4) 
(2.5) 



(2.6) 
(2.7) 



(2.8) 



Taking the ground state to have zero energy, the Hamiltonian operator in Eq. ( |2.5| ) (working in units where h = 1) 
can be split as H = Hq + Hi , where 



Ho = 2uj L \a 1 )(a 1 \ + 2uj L \a 2 ) {a 2 \ + u> L \b){b\ + u> d \d){d\ 
is approximately equal to the Hamiltonian of the molecular system, and 



(2.9) 



|oa»(6| 



— iUJLt 



+ H.c] 



H x = [n 6c |6)(c|e- <W£t + H.c] + [iUflai) 
+ [g(|ai) + |fl2»(c|e^ t +H.c] 

+(oji - 2oj L )\a 1 )(a 1 \ + (u> 2 - 2uj L )\a 2 )(a 2 \ + (uj b - uj L )\b)(b\ + (uj d - oj L )\d)(d\ 



(2.10) 



includes the interaction with the laser field plus corrections to H to reproduce the full molecular Hamiltonian. 

The first and second terms in Eq. (EETof) describe the interaction of the classical laser field with electric dipole 
moments of the one-photon transitions |c) — > \b) and \b) — > \ai), \a 2 ), respectively. The strengths with which these 
transitions are driven are characterized by the one-photon Rabi frequencies Of, c = \p\, c • El, and il ab — \p bai ■ El = 
\p ba2 ' E>l, where El is t he am plitude of the laser field. 

The third term in Eq. ( 2.10D describes the two-photon coupling of the laser field to the system with the two-photon 
Rabi frequency 



1 (^mcl^mai-^L 

2 lul — uj mc 



1 f-t"rncf^7na2-^ L 

2 LUl — UJjnc 



(2.11) 
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where El = \El\- This is due to transitions via the intermediate virtual levels labelled m here. 

Because of the external driving the elements of the system state matrix p satisfy equations of motion containing 
explicit time-dependent factors of the complex exponential type. These can be removed by moving to the interaction 
picture with respect to Hq. The remaining Hamiltonian Hi becomes 

Hi(t) = (wi 2 /2- A)|ai)(a 1 | 

+(-wi 3 /2- A)\a 2 ){a 2 \ 
+ (- A/2 -S)\b)(b\ 
+ [n ab (\ ai ) + \a 2 ))(b\+R.c.} 
+ [Q (M + |a 2 » <c| + H.c] 

+ [n 6c |6)(c|+H.c.]. (2.12) 

Although this is written as Hi(t) it is actually time- independent because of the judicious choice of Hq- Here A = 
2lul — u a is the detuning between the two-photon laser frequency 2ujl and the mean frequency of the upper levels 
relative to the ground level uj a = + lu 2 )/2. The one-photon detuning S — ljl — A/2 — u>t, — lu u /2 — u>b is the gap 
between the energy of level \b) and the half way position from the ground level |c) to the mean of the upper levels 
|eti) and \a 2 ). Moving to the interaction picture does not affect the irreversible terms so the new master equation is 

p = C ia p-i[Hi,p]. (2.13) 

The stationary solution satisfying p = can be found numerically and, in certain limits, analytically. We consider 
separately the case of two-photon coupling only, and one- and two-photon coupling. 



III. TWO-PHOTON COUPLING ONLY 



The case where the upper pair of levels is excited only by two -photon transitions via virtual intermediate levels is 
found by setting Q a b and Qt, c in the interaction Hamilltonian ( [2.12 ) equal to zero. The two-photon driving parametrized 
by Q is the only sort of driving mentioned in the experimental paper Q. 



1. Analytical Solution 



We first consider an analytical solution. This is possible in the weak-field limit where Q is much smaller than the 
decay rates in the system. For the experimentally relevant mutual polarizations p v = 1, p u — —1, the equations of 
motion are greatly simplified if we make the assumption that 7n — j v . That is, we assume that the decay rates of the 
upper levels on the ultraviolet and visible transitions are equal. We therefo re de fine a new parameter 7a — j u = j v . 

Under these assumptions, it is easy to show that the master equation (2.13) leads to the following steady-state 
values of the upper level populations and coherences 



PU 



p 22 



Q 2 



(A + Wl2 /2)2+ 7 2- 

Q2 

(A- wl2 /2)2+ 7 r 



Repi 



Q 2 [A 2 -(c 12 /2) 2 + 7q 2 ] 



[(A + c 12 /2) 2 + 7 2 ][(A-„ 12 /2) 2 + 7 2 ] 



(3.1) 
(3.2) 
(3.3) 



These are shown in Fig. ||(a) as a function of A. 

This analysis predicts that the populations and coherence exhibit peaks at A = ±u/i 2 /2, corresponding to the two- 
photon resonances of the laser field with the |c) — * |ai) and |c) — * \a 2 ) transitions. In Fig. ^(b), we plot the fluorescent 
intensity as a function of A for the p v — 1 and p u = — 1 transitions. It is seen that there are two peaks located at 
A = ±hu>i 2 , the amplitudes of which are not sensitive to p. The intensity is sensitive to p only about A = and 
can be almost completely suppressed for p v = 1 transitions. This confirms the earlier prediction by Agarwal jl(| that 
the two-photon excitation process involving the \a±) and \a 2 ) levels can lead to cancellation of spontaneous emission 
to the level \b). The cancellation of the fluorescence at A = also confirms the prediction by Berman JllJ that the 
suppression of the fluorescence can be explained in terms of dark states and coherent population trapping. 
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For p v = 1 the fluorescent intensity ( |2.3| ) can be written as 

I v = 2j v p ss , (3.4) 

where p ss — (s\p\s) is the population of the symmetric \s) = (|ai) + |a2) ) /V2 combination of the upper levels. 
The suppression of the fluorescence at A = indicates that the state |s) is almost unpopulated in the steady- 
state. This implies that the population is trapped between other molecular levels, including the antisymmetric state 
\a) = (\cii) — |a 2 )) /V%, with the state \s) being a dark state of the system. 



A. Numerical Results 



As noted before, in the experiment || three peaks were observed on the transitions with antiparallel dipole moments. 
However, as it is seen from Fig. ||(b), the weak-field theory does not predict three peaks for the p u = —1 transitions. 
The reason is that the magnitude of the coherence term p\2 is small compared to the magnitude of the population 
terms so that it is unable to build up a third peak in the middle. The coherence term is necessarily small because 
there is no detuning A at which both populations are large, and the coherence term is limited in magnitude by 

\P12\ 2 < /011/322- (3.5) 

To prove that the lack of a third peak is not a result of the assumptions made in deriving the analytical results we 



have also studied numerically the steady state of the master equation (2.13) with the one-photon Rabi frequencies set 



to zero. This can be done by calculating the equations of motion for the density matrix elements and using matrix 



inversion techniques. It can be done more easily using the direct symbolic representation of the master equation (2.13) 
which is possible in the quantum optics toolbox for matlab [Q. We find that, even in the strong field limit, and even 
with 7„ ^ 7„, it is not possible to produce a third peak in the fluorescence profile. 

From these analytical and numerical results we conclude that as well as the two-photon excitation process there 
must be some other processes involved in the dynamics of the system. The obvious candidate is a two-step one-photon 
process. 



IV. ONE- AND TWO-PHOTON COUPLING 



To include one-photon coupling we now consider the case where £lbc and Q a b are nonzero. To include two-photon 
coupling we actually do not need to have Q nonzero. That is because, as we will show, there is a regime in which 
level \b) acts as a virtual level with almost no real population. In this limit, the two step one-photon process becomes 
equivalent to a two-photon process. The relative strength of the two- and one-photon couplings is given by a parameter 
a — jb/ja, to be discussed later. Thus, for simplicity, we set Q = 0. 



A. Analytical Results 

To obtain analytical results we must consider the equations of motion for the density matrix elements. The 
master equation ( 2.13| ), in general, leads to a system of twenty five equations of motion for the density matrix ele- 



ments. Because of the assumption of large non-degeneracy between the intermediate levels \b) and \d), the coherences 
Pcd, Pdai , Pbd and pda 2 are n °t coupled to the driving field, and then the system of equations splits into two subsystems: 
one of seventeen equations of motion directly coupled to the driving field and the other of eight equations of motion 
not coupled to the driving field. It is not difficult to show that the steady-state solutions for the eight density matrix 
elements are zero and therefore we limit our considerations to the seventeen equations which, after applying the trace 
property (Trp =1), reduce to a system of sixteen coupled linear inhomogeneous equations. 



As in the case of Sec. Ill 1, for the physical parameters p v — 1, p u — — 1, the equations are simplified if 7„ = 7„. 
Under this assumption, and substituting 7 a for both j u and j v , the relevant density matrix elements obey the following 
coupled equations 

= ~2j a pn - itt ab {pbi - pu) , (4.1) 

P22 = -27a/0 2 2 - iV-ab {Pb2 - P2b) , (4.2) 
Pbb = -IbPbb + 7a (Pll + P22 + Pl2 + P2l) 
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- ifife c (p cb - Pbc) + iQab (Pbl - Plb) + i^ab (Pb2 - P2b) , 
Pec = IbPbb + 7d(l - Pll - P22 - Pbb - Pec) + i^bc (Pcb ~ Pbc) , 
Pl2 = - (27a + ^12) Pl2 - iQabPb2 + i^abPlb, 

Pb2 = - [(la + lb/2) -i(S— A/2 - W12/2)] p b2 - ifl b cPc2 - i^abPl2 + i^abPbb - i^abP22, 
Pbl = - [(la + lb/2) -i(S— A/2 + W12/2)] pbl - iflbcPcl - i&abPl2 + i^abPbb - iClabPu, 
Pbc = - [lb/ '2 -i(8+ A/2)] pbc - iftbc (Pec - Pbb) - iftab (Pic + P2c) , 
Pic = [-*(wi2/2 - A) - 7 a ] pic - iflabPbc - laPlc + i^bcPlb, 
P2c = [i(b)u/2 + A) - 7q] p2c + i^abPbc ~ laP2c ~ iftbcPlb- 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 
(4.10) 



To proceed further we make the weak-field assumption that f2 a 6 and f2 0C are small compared with the decay rates. We 
will show later that the same qualitative results can be obtained when this assumption, and the assumption j u = 7^, 
are relaxed. 

Under the weak field assumption we can order the matrix elements by how they scale with f2 ~ f^, f2 a b as shown 
in Fig. [|. The simplifications result from keeping only the lowest order terms in the above equations of motion for 
the state matrix elements. The steady-state solutions can then be obtained by setting the time derivatives to zero 
and solving the equations in the order as shown in Fig. til We find that the upper level populations are given by 



Pu 
P22 



ab be 



[(A/2 + ^+7 b 74][(A- Wl2 /2) 2 + 7a 2 ] 



fi 2 ,fi 2 

ab be 



[(A/2 + ( 5) 2 +7 2 /4][(A + a, 12 /2) 2 + 7a 2 ] 



(4.11) 
(4.12) 



This result predicts that, for large enough level splitting u>\2, the population of both of the upper pair states have 
two distinct peaks as a function of laser detuning A. This is illustrated in Fig. [|(a). The first peak is centered at 
A = ±0)12/2 for pu or P12 respectively. At this detuning the two-photon transition from |c) to \a\) or \a-i) is resonant, 
explaining the peak. The second peak is at A = —28. This is the resonance condition for the transition from |c) to 
I b), as seen in Fig. |l|. This central peak results from two stepwise one-photon transitions, the first populating level 
1 6) and the second exciting from \b) to \a). The populating of level \b) at this laser frequency is evident from the 
steady-state result 



Pbb 



fc"bc 



[(A/2 + 6^+ 7fc 2 /4] ■ 



(4.13) 



The upper-states coherence is considerably more complicated, and is given in full in the Appendix. From the 
denominators in the expression given there, it is evident that P12 may have many peaks and this is also illustrated 
in Fig. ^(a). To discover the physical meaning out of such a complicated expression, we consider the limit of large 
splitting where W12 is much larger than all other rates or frequencies. Wc then consider the behavior of pi2 at the 
positions of its peaks, and keep only the leading contributions there. It turns out that only one peak survives this 
simplification: 



Repi2 



ab be 



(u, 12 /2) 2 [(A/2 + 5f + ( lb /2f 



(4.14) 



From this expression it is evident that the coherence P12 also exhibits the resonance at A = —28, and its magnitude 
is comparable to the magnitude of population terms. This is possible because the populations pu and P22 both 
have peaks at A = —26, resulting from two-step one-photon transitions. Thus the inequality in Eq. ( |3.5| ) allows the 
coherence (4.14) to have a peak here also, unlike the case with only t wo-p hoton tran sitions. 

Assuming agai n tha t LO12 S> la->lb,8, the peaks in the populations ( 4.11 ) and ( 4.12 ) are well-separated Lorentzians. 
Then using Eq. (4.14), the fluorescent intensity for the ultraviolet and visible transitions can be approximated as 



167„ A Q 2 5 Q 2 C 



(l -Pu/v) 



(A-^ 12 /2) 2 + 7 2 2 (£ + A/2) 2 + ( 7fc /2) 2 (A + ^ 12 /2) 2 + 



(4.15) 



In this limit the fluorescent intensity contains three Lorentzians lo cated a t A = ±ui 2 /2 and A = —28. This is seen 
in the complete analytical solution for the fluorescence, from Eqs. (4.11), (4.12), and (5.1), plotted in Fig. Q(b). The 
amplitude of the peak at A = — 28 strongly depends on the mutual polarization of the dipole moments. The peak 
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is absent in the intensity observed in the visible region with p v = 1. For the fluorescent intensity observed in 
the uv region with p u = —1, the amplitude of the peak is enhanced. The strong dependence of the amplitude of 
the central peak on the mutual orientation of the molecular dipole moments is precisely the effect observed in the 
experiment. We emphasize again that the presence of the central peak in the fluorescent intensity results from the 
coupling of the driving laser to the one-photon transitions. 

B. Numerical results 

Having illustrated the role of the one- and two-photon excitations in the weak- field limit, we now find the fluorescent 
intensity without making any simplifying assumptions in our model. In this case it is not possible to obtain analytical 
solutions and therefore we use numerical methods to find stationary values of the density matrix elements of the 
system. Once again, this is easy using the symbolic representational power of the quantum optics toolbox for matlab 
. We first verify the correctness of the numerical technique by reproducing the weak-field analytical results. This 
is shown in Fig. §(c). 

In Fig. ||, we plot the fluorescent intensity for a strong driving field. It is seen that the fluorescent intensity exhibits 
the same behavior as that for the weak driving field, shown in Fig. |[ despite the fact that the solutions have been 
derived in different regimes. 

The experimentally observed fluorescent intensity was asymmetric about A = 0. There are few factors which 
could contribute towards the observed asymmetry. For example, the decay rates from the two upper levels to the 
intermediate levels cou ld be unequal. A simpler reason could be that the central peak is not exactly at A = 0. The 
analytical solution ( |2 .3| ) predicts the central peak to be at A = — 25 and the condition of 5 = implies that the energy 
of the level \b) is exactly half of the mean energy of the upper levels. There is no reason to expect this condition to 
be satisfied in the real molecule, and in fact it appears from the experimental results that 5 is positive. Fig. [7] shows 
the effect of a non-zero S on the fluorescence profile for a strong driving field. 

The relative magnitude of the central peak to the magnitude of the side peaks at A = ±wi2/2 depends on the ratio 
a = "fb/la- I n Fig- ^| we show the effect of a on the amplitude of the central peak in the fluorescent intensity on the uv 
transition. It is seen that the relative amplitude of the central peak increases with decreasing a (although the overall 
fluorescent intensity decreases) . The exact size of the central peak compared to the side peaks depends on f2 and W12 
as well as a. A small value of a, which results in a large central peak as observed in the experiment, is consistent with 
the fact that the decay rates of the intermediate levels are much smaller than the decay rates of the upper levels [^2) . 
When a increases the central peak becomes relatively smaller, and disappears completely for sufficiently large a. In 
this case the middle level is scarcely populated (because of its large decay rate) and the dynamics of the system are 
dominated by a two-photon process where the upper levels are directly populated from the ground level. Thus the 
large a limit is equivalent to considering only two-photon processes as in Sec. Ill, and it is not surprising that the 
spectrum contains only two peaks as found in that section. 

Finally, in Fig. ^| we show numerically that the results are not much affected if we relax our previous assumption 
that 7„ and 7„ are equal (with their value being denoted by j a )- For this plot we choose "f u and j v to be different by 
more than a factor of two. The numerical results in this figure, and all of the above figures, indicate that the existence 
of the third peak is a robust feature which does not depend upon fine tuning of the parameters in the model. 

V. SUMMARY 

We have modeled quantum interference effects in the intensity of the fluorescence emitted from a five-level molecular 
system, studied experimentally by Xia et al. J^] . We have presented an analytical solution for the fluorescent intensity, 
valid in the weak-field limit, and a numerical solution valid for arbitrary strengths of the driving field. We have been 
particularly interested in a theoretical explanation of the experimentally observed dependence of the number of peaks 
in the fluorescent intensity on the mutual orientation of the transition dipole moments. We have assumed that the 
molecular excitation is composed of a one-step, two-photon absorption process, and a two-step process involving 
the absorption of a single photon in each step. If the excitation is composed of only the two-photon processes, the 
fluorescent intensity consists of two peaks regardless the mutual orientation of the molecular dipole moments. With 
the two-step, one-photon processes included, the intensity consists of two peaks on transitions with parallel dipole 
moments and three peaks on transitions with antiparallel dipole moments. This latter case is in excellent agreement 
with the experimental observation The variation of the number of peaks with the mutual polarization of the 
dipole moments is a very clear demonstration of quantum interference in spontaneous emission. 
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APPENDIX 

The complete analytical solution to the upper level coherence in the weak driving limit with 7„ = 7„ = j a is 

R0P12 = (tt ab n bc ) 2 [ lb - /a (5 - A/2 + wi 2 /2)(A - wia/2) + ibla 2 {la + 7b)/2 

+2 7a 2 (<5 - A/2 + w 12 /2)(A/2 + <J) - 7 a(2 7 a + 7b) (A - w 12 /2)(A/2 + S)/2 
+uj 12 (5 - A/2 + wi 2 /2)(A - wi 2 /2)(A/2 + <J) + ^i 27a (2 7a + 7b)(A/2 + 5)/2 
-^I27a7b(<5 - A/2 + wi 2 /2)/2 + wi 27b (2 7a + 7b )(A - wi 2 /2)/4] 

- { [(<5 - A/2 + Wl2 /2) 2 + (2 7a + 7b ) 2 /4] [(A - uj 12 /2) 2 + 7a 2 ] 

[( ( 5 + A/2) 2 + 7 2 /4][47 a 2 +^i 2 2 ]} 

+ (n ab n bc ) 2 [ 7 i,7 a (<5 - A/2 - wi 2 /2)(A + wi 2 /2)/2 + 7 b7a 2 (2 7a + 7b)/2 
+2 7a 2 (<5 - A/2 - wi 2 /2)(A/2 + <$) - 7a (2 7a + 7b 2)(A + w 12 /2)(A/2 + 5)/2 
+w 12 (<5 - A/2 - wi 2 /2)(A + wi 2 /2)(A/2 + (5) + ^ 12 (2 7a + 7 b) 7 a(A/2 + 5)/2 
-^I27a7b('5 - A/2 - wi 2 /2)/2 + ^i 27b (2 7a + 7b )(A + wi 2 /2)/4] 

- { [(J - A/2 - Wl2 /2) 2 + (2 7a + 7b ) 2 /4] [(A + Wl2 /2) 2 + 7a 2 ] 
[(<5 + A/2) 2 + 7b 2 /4][4 7a 2 + Wl2 2 ]} 

, ^M-la^la + 7b) -(6- A/2 + uWgMg] 

+ [(5 - A/2 + Wl2 /2)2 + (2 7a + 76)2/4] [((5 + A/2) 2 + 7 2 /4] [4 7a 2 + Wl2 2] 

, »q b »L[-7a(27a + 7b) - (6 - A/2 - aWgVigj 

+ [(6 - A/2 - wi 2 /2)2 + (2 7a + 76)2/4] [(5 + A/2)2 + 7 2 /4] [4 7a 2 + Wl2 2 ] ' 
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FIG. 1. Energy level structure and couplings of the molecular system. 
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FIG. 2. Analytical results for 2-photon coupling only, (a) shows populations (pn peaked to the left, p22 peaked to the right) 
and coherence (5ftpi2 , below the axis) of the excited states, (b) shows the total fluorescent intensity in units of the spontaneous 
emission rate. The solid line shows the intensity on the ultraviolet transition (p u = — 1), the dashed line the intensity on the 
visible transition (p v = 1), and the dash-dot line the hypothetical intensity for a transition with orthogonal dipole moments 
(p — 0). The parameters are Q, a b = ilbc = 0, Q = 10~ 4 ,tJi2 = 6, 8 — 0, ~f u — "fv = 0.5,76 = 1. The two-photon detuning A is 
plotted in units of j u + y v . 
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FIG. 3. Diagram showing the method for solving the steady state Master equation under the weak field assumption. The 
symbols on the right represent the order of the matrix elements. 
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FIG. 4. Weak field results for one- and two-photon driving, (a) shows analytical populations and coherences of the 
the excited states, as in Fig. ^(a). (b) shows the analytical results and (c) the numerical results for the total fluores- 
cent intensity in units of the spontaneous emission rate. The different line styles are as in Fig. ^|(b). The parameters are 
tlab = Qbc = 0.01, Q = 0, W12 = 6, 8 = 0, -y u = 7„ = 0.5, 7b = 1. The two-photon detuning A is plotted in units of j u + j v . 




FIG. 5. Total fluorescent intensity in units of the spontaneous emission rate for strong one- and two-photon coupling. The 
different line styles are as in Fig. |^(b). The parameters are Q a b = Qbc = 1, Q = 0, W12 = 6, 8 = 0, 7 U = j v = 0.5, 75 = 0.15. The 
two-photon detuning A is plotted in units of j u + "fv ■ 




FIG. 6. Total fluorescent intensity in units of the spontaneous emission rate for strong one- and two-photon coupling, 
with non-zero detuning S of the intermediate level |6). The different line styles are as in Fig. |^(b). The parameters are 
Q a b = fibc = 1, Q = 0, oji2 = 6,5 = 0.3, 7u = 7t, = 0.5, fb = 0.15. The two-photon detuning A is plotted in units of 7^ + j v . 
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FIG. 7. Total fluorescent intensity in units of the spontaneous emission rate for strong one- and two-photon coupling. Only 
the ultraviolet (p u = — 1) transition is plotted, but 71, = wy u is varied. The values of a for the three curves are, from top to 
bottom, 2, 0.3, and 0.02. The other parameters are £l a b = f2tc = 1, Q = 0,cji2 = 6, S = 0, j u — "fv = 0.5. The two-photon 
detuning A is plotted in units of 7^ + 7^. 




FIG. 8. Total fluorescent intensity in units of the spontaneous emission rate for strong one- and two-photon coupling, with 
non-equal decay rates on the ultraviolet and visibile transitions. The different line styles are as in Fig. ^|(b). The parameters 
are Q a b ~ Qbc = l,Q — 0,a>i2 = 6,5 = 0,7 U = 0.7, 7„ = 0.3,7b = 0.15. The two-photon detuning A is plotted in units of 
7« + lv ■ 
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